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ABSTRACT 

We study the hard-X-ray spectral properties of ten highly luminous radio-quiet (RQ) active galactic nuclei 
(AGNs) at z = 1.3 — 3.2, including new XMM-Newton observations of four of these sources. We find a sig- 
nificant correlation between the normalized accretion rate (L/Lndd) and the hard-X-ray photon index (T) for 
35 moderate-high luminosity RQ AGNs including our ten highly luminous sources. Within the limits of our 
sample, we show that a measurement of T and Lx can provide an estimate of L/^Edd an d black-hole (BH) mass 
(Mbh) with a mean uncertainty of a factor of ^ 3 on the predicted values of these properties. This may provide 
a useful probe for tracing the history of BH growth in the Universe, utilizing samples of X-ray-selected AGNs 
for which L/Lgdd an d Mbh have not yet been determined systematically. It may prove to be a useful way to 
probe BH growth in distant Compton-thin type 2 AGNs. We also find that the optical-X-ray spectral slope 
(a ox ) depends primarily on optical-UV luminosity rather than on L/Lgdd in a sample of RQ AGNs spanning 
five orders of magnitude in luminosity and over two orders of magnitude in L/Lgdd- We detect a significant 
Compton-reflection continuum in two of our highly luminous sources, and in the stacked X-ray spectrum of 
seven other sources with similar luminosities, we obtain a mean relative Compton reflection of R = 0.9^q 5 
and an upper limit on the rest- frame equivalent width of a neutral Fe Ka line of 105 eV. We do not detect 
a significant steepening of the X-ray power-law spectrum below rest-frame 2 keV in any of our highly lumi- 
nous sources, suggesting that a soft-excess feature, commonly observed in local AGNs, either does not depend 
strongly on L/L^dd, or is not accessible at high redshifts using current X-ray detectors. None of our highly 
luminous sources displays X-ray flux variations on timescales of ~ 1 hr, supporting the idea that the timescale 
of X-ray variability depends inversely on Mbh an d does not depend on L/Ledd- 
Subject headings: galaxies: active - galaxies: nuclei - X-rays: galaxies - quasars: emission lines 



1. INTRODUCTION 

The X-ray spectrum of an unobscured, radio-quiet (RQ) 
AGN in the ^2-100 keV energy band is best characterized 
by a single power-law continuum of the form N(E) °= E~ r , 
where T, hereafter the photon index in the ^2-100 keV en- 
ergy band, typically lies in the range ^1.5-2.5. A corona of 
hot electrons is assumed to produce the hard-X-ray emission 
via Compton upscattering of UV-soft-X-ray photons from the 
accretion disk, and F is predicted to be only weakly sensitive 
to large changes in the electron temperature and the optical 
depth in the corona (e.g., Haardt & Maraschi 1991; Zdziarski 
et al. 2000; Kawaguchi et al. 2001). The relatively narrow 
range of F values in RQ AGNs has been reported in numerous 
studies (e.g., Nandra & Pounds 1994, Reeves & Turner 2000; 
Page et al. 2005; Shemmer et al. 2005; Vignali et al. 2005; 
Just et al. 2007), and typically no strong dependence of F on 
redshift or luminosity has been detected across the widest pos- 
sible ranges of these parameters. On the other hand, a strong 
anticorrelation between F and the full width at half-maximum 
intensity (FWHM) of the broad emission-line region (BELR) 
component of Hj3 has been found, first by Brandt et al. (1997). 

The remarkable dependence between X-ray and optical 
spectroscopic properties has been suggested to arise from a 
more fundamental correlation between F and the accretion 
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rate (e.g., Brandt & Boiler 1998; Laor 2000). A high accre- 
tion rate is expected to soften (steepen) the hard-X-ray spec- 
trum by increasing the level of disk emission, resulting in 
the production of softer photons, which increase the Comp- 
ton cooling of the corona. Using recent scaling relations for 
the BELR size, luminosity, and the width of the broad Hj3 
emission line from reverberation-mapping studies, it is clear 
that the normalized accretion rate (i.e., L D0 i/£Edd, nerea ft er 
L/LEdd, where Lboi is the bolometric luminosity) is propor- 
tional to FWHM(Hj3)~ 2 , at least for low-moderate luminos- 
ity AGNs in the local universe (e.g., Kaspi et al. 2000). 

Subsequent X-ray studies of nearby (z 0.5) unob- 
scured RQ AGNs have confirmed the Brandt et al. (1997) 
r-FWHM(Hj3) anticorrelation (e.g., Leighly 1999; Reeves 
& Turner 2000; Porquet et al. 2004; Piconcelli et al. 2005; 
Brocksopp et al. 2006), and others have found significant 
correlations between F and L/Ledd (e.g., Lit & Yu 1999; 
Porquet et al. 2004; Wang et al. 2004; Bian 2005). How- 
ever, all these studies were not able to disentangle the strong 
FWHM(Hj3)-L/LEdd dependence. Recently, Shemmer et al. 
(2006; hereafter S06) have suggested that this degeneracy can 
be removed if highly luminous sources are included in the 
analysis. This can be achieved by obtaining high-quality near- 
IR spectroscopy of the H/3 spectral region, to obtain L/Lndd 
(e.g., Shemmer et al. 2004), as well as accurate measurements 
of r using XMM-Newton and Chandra for highly luminous 
AGNs found at 1 < z < 3. 

Within the limits of their sample of 30 sources, spanning 
three orders of magnitude in luminosity, S06 have shown 
that r does not depend on luminosity or black-hole (BH) 
mass (Mbh)- They have also shown that the F values of the 



2 



SHEMMER ET AL. 



five highly luminous sources in their sample, while consis- 
tent with the values expected from their normalized accretion 
rates (L/LgddX are significantly higher than expected from the 
widths of their broad H/3 emission lines. This has enabled, 
for the first time, breaking of the FWHM(H/3)-L/L E dd degen- 
eracy and has provided evidence that F depends primarily on 
the accretion rate. However, the number of highly luminous 
sources was small, and as explained below, some uncertainties 
remained. 

In this work, we double the number of highly luminous 
sources at high redshift and reinforce the S06 results. We 
show that r can be considered a reliable accretion-rate indi- 
cator for moderate-high luminosity RQ AGNs, and that the 
combination of F and X-ray luminosity may provide a useful 
probe for tracing the history of BH growth in the Universe. 
We also discuss the X-ray spectral and temporal properties of 
luminous, high-accretion rate AGNs at high redshift as well 
as the dependence of the optical-X-ray spectral energy distri- 
bution (SED) on luminosity and L/Lndd- In §|2]we describe 
our sample selection, and present the new observations and 
their analysis. Our results are presented and discussed in §[3] 
where we focus on the correlation between the hard-X-ray 
photon index and the normalized accretion rate in RQ AGNs 
and its implications for probing BH growth in the Universe. 
A summary of our main findings is given in §|4] Throughout 
this work we consider only RQ AGNs to avoid any contribu- 
tion from jet-related emission to the X-ray spectra. Luminos- 
ity distances are computed using the standard cosmological 
model with D. A = 0.7, D. M = 0.3, and// = 70 km s _1 Mpc~ ! . 

2. SAMPLE SELECTION, OBSERVATIONS, AND DATA ANALYSIS 

2.1. Sample Selection 

Our sample is composed of the 30 sources stud- 
ied in S06 as well as five highly luminous sources 
(46 < log[vL v (5100 A)] < 48) whose hard-X -ray spectra are 
analyzed in this work for the first time. Motivated by the S06 
hypothesis that only highly luminous AGNs may break the 
FWHM(Hj3)-L/LEdd degeneracy, we have primarily searched 
for highly luminous AGNs that have high-quality Hj3 and 
hard-X-ray spectroscopy. Our search of the literature and 
archive yielded ten such sources, five which were studied in 
S06 and five that are analyzed below. Our new core sample 
of 35 sources includes unabsorbed, type 1 RQ AGNs with 
44 < log[vL v (5100 A)] < 48 (i.e., a moderate-high luminos- 
ity range; see Fig. [T). Rest-frame optical data for this sam- 
ple was obtained from Neugebauer et al. (1987), Boroson & 
Green (1992), Nishihara et al. (1997), Mcintosh et al. (1999), 
Shemmer et al. (2004), and Sulentic et al. (2006). Hard- 
X-ray data were obtained from Page et al. (2004b), Picon- 
celli et al. (2005), S06, and from this work (see below). All of 
our sources were selected to have high-quality XMM-Newton 
hard-X-ray spectra (except for HE 0926—0201 that was 
serendipitously observed with Chandra). We have intention- 
ally limited our sample to the moderate-high luminosity range 
in order to minimize potential effects of optical and X-ray 
variability (since the optical and X-ray data are not contem- 
poraneous) as well as potential spectral complexities due to 
Compton reflection, as discussed further below (see also S06 
for more details). In light of our selection criteria, we caution 
that our sample is neither complete nor fully representative of 
the AGN population as a whole; our main results may there- 
fore be subject to selection biases over the wide range of AGN 
properties, which were not fully explored in this work. 
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FIG. 1. — Luminosity histogram of our core sample of 35 AGNs. The 
shaded region marks the ten highly luminous sources at z = 1.3 — 3.2 (the 25 
moderate-luminosity sources are at z < 0.5). 

2.2. New XMM-Newton Observations 

We have performed X-ray spectral imaging observations of 
four new sources from the Shemmer et al. (2004) sample with 
XMM-Newton (Jansen et al. 2001); a log of these observa- 
tions is presented in Table[T] These sources were selected 
for XMM-Newton observations for being luminous, high- 
accretion rate sources, predicted to have high X-ray fluxes, 
and having relatively low Galactic column densities. The 
data were processed using standard XMM-Newton Science 
Analysis System 5 v6.5.0 tasks. The event files of all the ob- 
servations were filtered to remove periods of flaring activ- 
ity in which the count rates of each MOS (pn) detector ex- 
ceeded 0.35 (1.0) counts s _1 for events having E > lOkeV. 
The time lost due to flaring in each observation varied be- 
tween l%-70% of the entire observing time; the net expo- 
sure times in TableQ] reflect the filtered data. The X-ray 
spectra of the quasars were extracted from the images of all 
three European Photon Imaging Camera (EPIC) detectors us- 
ing apertures with radii of 30". Local background regions 
were at least as large as the source regions. The spectra were 
grouped with a minimum of 20 counts per bin, except for the 
spectrum of Q 1318 — 113 which was grouped with a mini- 
mum of 50 counts per bin. Joint spectral fitting of the data 
from all three EPIC detectors for each source was performed 
with XSPEC vl 1.3.2 (Arnaud 1996). We employed Galactic- 
absorbed power-law models at rest-frame energies ^ 2 keV, 
corresponding to ^ 0.5-0.6 keV in the observed frame of the 
sources, where the underlying power-law hard-X-ray spec- 
trum is less prone to contamination due to any potential soft 
excess emission or absorption. In each fit, the photon indices 
in the spectra of all three EPIC detectors were tied to a single 
value, while the normalizations were free to vary. The best-fit 
r values, power-law normalizations, and % 2 values from these 
fits are given in columns (2), (3), and (4) of Table|2] respec- 
tively, and the data, their joint, best-fit spectra, and residuals 
appear in Fig. [2] 

We also searched for intrinsic absorption in each source by 
jointly fitting the spectra with a Galactic-absorbed power law 
model including an intrinsic (redshifted) neutral-absorption 
component with solar abundances in the same energy range 
quoted above. No significant intrinsic absorption was de- 
tected in any of the sources; upper limits on intrinsic Nn val- 
ues appear in column (5) of Table[2] Each panel of Fig.|2]in- 
cludes a F-Nr confidence-contour plot from this fitting for 
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TABLE 1 
XMM-Newton OBSERVATION LOG 













Observation 


Net Exposure Time (ks) / Source Counts 


AGN 


RA (J2000.0) 


DEC (J2000.0) 


z a 


Nn b 


Start Date 


MOS1 MOS2 


PN 


LBQS 0109+0213 


01 12 16.91 


+02 29 47.6 


2.349 


3.25 


2007 Jan 8 


33.7 /211 33.8/218 


26.3 /554 


2QZJ023805. 8-274337 


02 38 05.80 


-27 43 37.0 


2.471 


1.73 


2006 Dec 22 


34.6/206 34.0/183 


23.9 / 586 


Q 1318-113 


13 21 09.38 


-11 39 32.0 


2.306 


2.80 


2006 Dec 28 


47.7 / 853 46.9 / 876 


33.2/2621 


SBS 1425+606 


14 26 56.18 


+60 25 50.9 


3.202 


1.58 


2006 Nov 12 


20.6 / 232 20.9 / 255 


7.0 / 304 



NOTE. — The net exposure time represents the L1VET1ME following the removal of flaring periods. The net source counts are in the 0.2-10.0 keV band. 
1 Systemic redshift measured from the optical emission lines and obtained from Shemmer et al. (2004). 
' Neutral Galactic absorption column density in units of 10 20 cm" 2 obtained from Dickey & Lockman (1990). 



each source. By applying F-tests between the models in- 
cluding intrinsic absorption and those that exclude it, we 
found that none of the spectra require an intrinsic absorp- 
tion component. The remarkably flat hard-X-ray spectrum of 
LBQS 0109+0213 (with T = 1.23; see Table EJ motivated 
an alternative modeling, searching for an indication of partial 
covering. We found that when the spectrum is fitted with a 
redshifted partial-covering fraction absorber (ZPCFABS model 
in XSPEC) and a Galactic-absorbed power law, the covering 
fraction is consistent with zero and the fit is not improved with 
respect to a Galactic-absorbed power law model. Therefore, 
the unusually flat hard-X-ray spectrum of LBQS 0109+0213 
cannot be explained as being due to partial covering. 

2.3. The Chandra Spectrum of HE 0926-0201 

HE 0926-0201 is a highly luminous RQ AGN at z = 1 .682 
for which near-infrared spectroscopy of the Hj3 region is pre- 
sented in Sulentic et al. (2006). This source is serendipi- 
tously detected in an 8.8 ks Chandra/ ACIS-S observation from 
2002 February 20 (ID 3139) with 126 source counts in the 
0.5-8 keV observed-frame band. We analyzed the Chandra 
observation using standard CIAO 6 V3.2 routines. The spec- 
trum was extracted using PSEXTRACT from a circular aper- 
ture with a radius of eight pixels, corresponding to ~4", and 
the events were grouped to have a minimum of 10 counts 
per bin (the extremely low Chandra background was deter- 
mined from a source-free annular region with inner and outer 
radii of 10" and 25", respectively, centered on the source). 
The spectrum was fitted over the £ 2 keV rest-frame energy 
band using XSPEC with a model including a power law and 
a Galactic-absorption component with a column density of 
A/h=3.17 x 10 20 cm -2 . The data were also fitted with an ad- 
ditional intrinsic-absorption component, similar to the proce- 
dure described above for the XMM-Newton observations, and 
we found that the data do not warrant any additional absorp- 
tion. The best-fit spectral parameters as well as the upper limit 
on the intrinsic absorption are included in Table [2] and the 
Chandra spectrum and best-fit model appear in Fig. [3] 

3. RESULTS AND DISCUSSION 
3.1. Spectral and Temporal Properties 
3.1.1. Optical Luminosities and FWHM( Hp ) 

Basic optical spectroscopic properties of the new high- 
redshift sample are given in Table[2] The monochromatic lu- 
minosity at a rest-frame wavelength of 5100 A [vL v (5100 A)] 
is given in column (6), and FWHM(H/3) is given in column 
(7); except for HE0926-0201, for which the optical data 
were obtained from Sulentic et al. (2006), these data were 
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obtained from Netzer et al. (2007), who recently presented 
new and improved spectroscopic measurements for all the 
Shemmer et al. (2004) sources. The Mbh and L/Lgdd val- 
ues in columns (8) and (9), respectively, were determined as 
in S06, using the vL v (5100A) and FWHM(H/3) values given 
in columns (6) and (7), respectively, and based on the recent 
reverberation-mapping results of Peterson et al. (2004) and 
the Kaspi et al. (2005; hereafter K05) BELR size-luminosity 
relation (see also Kaspi et al. 2000); the general expression we 
use for L/L Edd is of the form L/L Edd °c L°- 3 FWHM(Hj3)- 2 . 
We note that the K05 relation relies on a sample of AGNs 
having luminosities up to vL v (5100A) «10 46 ergs s _1 , and 
extrapolating it to higher luminosities, such as those of ten 
of the sources presented here, is somewhat uncertain; a 
reverberation-mapping effort is underway to test the valid- 
ity of such extrapolations (see e.g., Kaspi et al. 2007). In 
addition, Bentz et al. (2006, 2007) have recently suggested 
that subtraction of host-galaxy starlight from the AGN op- 
tical continuum may result in a somewhat flatter slope (e.g., 
a = 0.54 compared with a = 0.69 in K05) for the BELR size- 
luminosity relation across the entire K05 luminosity range. 
We test the possibility of using a flatter slope below, but 
we note that the Bentz et al. measurements were performed 
for sources with vL v (5100A) < lO^ergss" 1 , while all the 
sources studied here have vL v (5100A) ^ 10 44 ergss _1 (with 
a median luminosity of vL v (5100A) ~10 45 ergs s , Fig.[T|). 
We also expect a much smaller fractional amount of host con- 
tamination in the most luminous sources in our sample. In 
light of all this, we have computed Mbh values using the K05 
slope across our entire luminosity range. Moreover, since in 
this work we employ nonparametric statistical ranking tests, 
our main results are not significantly sensitive to the precise 
value of the slope of the BELR size-luminosity relation in the 
range a — 0.54 — 0.69. To derive the normalized accretion 
rates, we have employed the luminosity-dependent bolomet- 
ric correction method of Marconi et al. (2004); the bolometric 
correction factor for the monochromatic luminosity at 5 100 A 
[i.e., L b oi/vL v (5100A)] is —6-8 (~ 5) for the moderate- 
luminosity z < 0.5 (high-luminosity z= 1 .3 — 3.2) sources in 
this work (see S06 for more details). 

3.1.2. Optical-X-ray Spectral Slopes 

The optical-X-ray spectral slopes (a ox ) in column (10) of 
Table [2] are defined as 

«ox= log(/2kev//2500A)/ lo g( V 2kev/v 25 ooA)' where /2keV 

and /2500 A ^ tne f ux densities at 2 keV and 2500 A, respec- 
tively. The 0£ O x values were derived using the photon indices 
and fluxes in columns (2) and (3), respectively, and the optical 
luminosities in column (6), assuming a UV continuum of the 
form f v oc v~° 5 (Vanden Berk et al. 2001). The Oo X values 
for all of our sources are consistent with the expected values, 
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FIG. 2. — Data, best-fit spectra, and residuals of our new XMM-Newton observations. Open circles, filled squares, and open squares represent the EPIC pn, 
MOS 1 , and MOS2 data, respectively. Solid lines represent the best-fit model for each spectrum, and the thick line marks the best-fit model for the pn data. The 
data were fitted with a Galactic-absorbed power-law model above a rest-frame energy of ~2keV, and then extrapolated to 0.3 keV in the observed frame. The 
X residuals are in units of o~ with error bars of size 1. The insets show 68%, 90%, and 99% confidence contours for T and Nu, when the data are fitted with an 
additional neutral intrinsic-absorption component. 

fleeted from neutral material was carried out by fitting all 
our XMM-Newton spectra in the ^ 2 keV rest-frame energy 
range with XSPEC, employing a Galactic-absorbed power- 
law and a Compton-reflection continuum model (i.e., the 
PEXRAV model in XSPEC; Magdziarz & Zdziarski 1995) si- 
multaneously with a redshifted Gaussian emission line model 
(using the ZGAUSS model in XSPEC); the Gaussian rest- 
frame energy and width were fixed at E = 6.4 keV and 
a = 0.1 keV, respectively. We also included in the analy- 
sis the five luminous, high-redshift sources from S06 that 
have XMM-Newton spectra and H/3 spectroscopy, namely 
PG 1247+267, Q 1346-036, PG 1630+377, PG 1634+706, 
and HE 22 17— 2818; the luminosities and accretion rates of 
these sources are comparable to those of the sources pre- 



given their optical luminosities (e.g., Steffen et al. 2006). 

3.1.3. Compton Reflection and Fe Ka Emission 

The relatively high redshifts of our sources allowed us 
to search for a Compton-reflection continuum as well as 
Fe Ka emission in their XMM-Newton spectra. A Compton- 
reflection continuum may be observed in AGNs within the 
rest-frame ~ 7 — 60keV energy range, peaking at rest-frame 
~30 keV, and it presumably originates from reflection of hard- 
X-ray photons off the relatively colder outskirts of the ac- 
cretion disk and/or the torus. We also searched for signa- 
tures of a neutral narrow Fe Ka emission line at rest-frame 
6.4 keV, as this is expected to appear in conjunction with a 
Compton-reflection continuum. The search for continua re- 
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FIG. 3.— The Chandra spectrum of HE 0926-0201. The solid line rep- 
resents the best fit Galactic-absorbed power law model. The % residuals are 
in units of o~ with error bars of size 1. The inset shows 68%, 90%, and 99% 
confidence contours for T and N}i, when the data are fitted with an additional 
neutral intrinsic-absorption component. 

sented in this work. Table[3]lists the best-fit parameters from 
these fits. Column (3) gives the rest-frame equivalent width 
(EW; throughout the paper, EWs refer to rest-frame values) 
of the Fe Ka emission line and column (4) gives the relative- 
reflection component (R) of the Compton-reflection contin- 
uum expressed as R = £l/2n, where £2 is the solid angle sub- 
tended by the continuum source. 

Except for the cases of PG 1247+267 and PG 1634+706, 
F-tests carried out between the results of the 
PEXRAV+ ZGA USS fits and the Galactic-absorbed power- 
law fits (§ !2.21 i have indicated that the new (reflection) model 
did not improve the fits with respect to the corresponding 
Galactic-absorbed power-law fits. Table [3] shows that, except 
for PG 1630+377, we have not detected any neutral Fe Ka 
emission in any of the sources, and instead, we have placed 
upper limits on the EWs of such emission. Although the 
EW(Fe Ka) for PG 1630+377 has a nonzero lower limit, and 
the EW value we found is consistent with the one found for 
this source by Jimenez-Bailon et al. (2005), we argue that 
the XMM-Newton data for PG 1630+377 do not warrant 
the existence of a reflection spectrum based on the F-test 
described above. 

We have also tested whether the XMM-Newton spectra of 
our sources exhibit hard-X-ray emission reflected from ion- 
ized material. This was performed by employing the same 
XSPEC fitting routine described above for the case of reflec- 
tion from neutral material, except that the PEXRAV model 
was replaced with the PEXRIV model (Magdziarz & Zdziarski 
1995), and the energy of the Fe Ka line was free to vary in 
the 6.7-6.97 rest-frame energy range to account for an ion- 
ized iron line. The results of these fits, i.e., R parameters and 
EW(Fe Ka) values, are consistent with those presented in Ta- 
ble [3] for the case of reflection from neutral material (even 
though the energies of the Fe Ka lines were required to be 
somewhat higher). 

The non-detections of Fe Ka emission lines in the highly 
luminous sources under study are in accord with previous 
studies claiming that the strength of such emission lines de- 
creases with increasing luminosity (aka the 'X-ray Baldwin 
Effect'; e.g., Iwasawa & Taniguchi 1993; Page et al. 2004a; 
Zhou & Wang 2005; Bianchi et al. 2007). However, the rela- 
tively high upper limits we obtain for our EW(Fe Ka) values 
cannot rule out the possibility that such an 'X-ray Baldwin Ef- 
fect' is very weak or does not even exist (e.g., Jimenez-Bailon 
et al. 2005; Jiang et al. 2006). Given the high luminosities 
of our sources, the expected EW(Fe Ka) values from such 
an EW-L relation are of order — 50eV, while in most cases, 
our upper limits are considerably higher (see Tableland e.g., 
Bianchi et al. 2007). 

We detect statistically significant Compton-reflection con- 



tinua only in PG 1247+267 and PG 1634+706 (Table 0. 
The R value we find for the first of these sources is con- 
sistent with the one found by Page et al. (2004b) for the 
same XMM-Newton spectrum. However, in contrast with 
our results, Page et al. (2004b) have not detected a significant 
Compton-reflection feature in the same XMM-Newton spec- 
trum of PG 1634+706 that we have analyzed, while the upper 
limit we obtain on the EW(Fe Ka) for the source is consistent 
with their finding; we also note that Nandra et al. (1995) have 
not detected any reflection features in an ASCA spectrum of 
the source. In addition, Page et al. (2004b) have detected an 
Fe Ka line in PG 1247+267 while we have not. Interestingly, 
inspection of Fig. 1 of Page et al. (2004b) suggests the exis- 
tence of a reflection component in both sources. An alter- 
native spectral fitting, including a thermal component and a 
power-law spectrum for each source, was recently employed 
by Ruiz et al. (2007); their study finds a soft excess and no 
reflection component in both sources. These partial discrep- 
ancies may be a consequence of the different modelings and 
different energy ranges used in the different studies. The ap- 
parently non-physical R values we obtain for PG 1247+267 
and PG 1634+706 (i.e., £2 > 4n steradians) may be due to the 
effects of general relativistic light bending (e.g., Fabian et al. 
2002; Fabian & Vaughan 2003). 

To improve our ability to detect or constrain the strengths 
of any reflection features, we jointly fitted the XMM-Newton 
spectra of our nine high-redshift sources from Table [3] In 
the joint fitting, we used the PEXRAV+ZGAUSS models as be- 
fore and tied all the R values and normalizations to a single 
value (the flux level of each source was controlled by assign- 
ing to it a scaling factor that was allowed to vary freely, and 
the photon index for each source was also free to vary). We 
ran the joint fitting process four times. In the first run, all 
nine sources were considered. In the second run, we removed 
the two sources from Table[3]in which a significant Compton- 
reflection continuum was detected. In the third run, we re- 
moved all sources in which a non-zero R value was detected, 
and in the fourth run we removed PG 1630+377 in which an 
Fe Ka line was detected. 

The best-fit values of EW(Fe Ka) and R from the joint- 
fitting process appear in Table @] where we also provide the 
total number of counts used in the fit and the mean redshift of 
the sources considered in each run. By comparing the results 
in Table @] with those in Table [5J one can see that when all 
sources are considered (run I), a large value of R is detected, 
while a tight constraint is placed on the mean EW(Fe Ka); 
this result is dominated by the properties of PG 1247+267 and 
PG 1634+706 that comprise -60% of the counts. In run II, 
one can see that despite the removal of all sources in which 
a Compton-reflection continuum was detected, a mean reflec- 
tion component is detected, although with a relatively small 
value of R. In run III, the Fe Ka detection is probably due 
to the fact that PG 1630+377 comprises -30% of the total 
counts; no significant Compton reflection is detected in this 
case. In the final run, no significant mean reflection emission 
is detected; the constraint on R is relatively tight while the 
constraint on EW(Fe Ka) is rather weak. 

These results do not provide a clear picture for the depen- 
dence of the reflection spectrum on luminosity or the accre- 
tion rate. Most of our high-redshift sources do not exhibit 
significant reflection components, as might be expected given 
their very high luminosities. Only two of these sources show 
a significant Compton 'hump', and a strong Fe Ka line is de- 
tected in another. In addition, significant detections of Comp- 
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TABLE 2 

Best-Fit X-Ray Spectral Parameters and Optical Properties 



logvL v (5100A) FWHM(Hj8) logM BH 



AGN 
(1) 


r 

(2) 


f v (l keV) a 
(3) 


X 2 /v 
(4) 


N H b 
(5) 


(ergss ') 
(6) 


(kins" 1 ) 
(7) 


(Me) 
(8) 


i/^Edd 

(9) 


(10) 


LBQS 0109+0213 
2QZJ023805. 8-274337 
Q 1318—113 
SBS 1425+606 
HE 0926-0201 


1 7^+0.12 
1 - ZJ -0.12 
9 13+0.16 
Z - 1J -0.15 

1 96+° 07 
1 76+ 014 

9 9O+0.36 


6 1 +0 ' 8 
oa -Q.7 

7 8+0- 8 
'■ S -0.8 

24.8+J2 


70/53 
47/39 
70/62 
38/39 
2/10 


< 1.46 

< 3.71 
<0.38 

< 2.33 

< 1.72 


46.8 
46.6 
46.9 
47.4 
47.0 C 


7959 
3403 
4665 
4964 
5100 c 


10.4 
9.5 
10.0 
10.4 
10.1 


0.1 
0.5 
0.3 
0.4 
0.3 


-1.89 
-1.66 
-1.64 
-1.82 
-1.80 



NOTE. — The best-fit photon index, normalization, and % 2 were obtained from a Galactic-absorbed power-law model. Errors represent 90% confidence 
limits, taking one parameter of interest (A% — 2.71). The optical data in columns 6 and 7 were obtained from Shemmer et al. (2004) and Netzer et al. (2007) 
except for HE 0926-020 1 . 

a Power-law normalization given as the flux density at an observed-frame energy of 1 keV with units of 10~ 32 ergs cm -2 s -1 Hz -1 ; except for HE 0926— 0201, 
this refers to the pn data, taken from joint fitting of all three EPIC detectors with the Galactic-absorbed power-law model. 

Intrinsic column density in units of 10 22 cm -2 . Upper limits were computed with the intrinsically absorbed power-law model with Galactic absorption, and 
represent 90% confidence limits for each value. 
c Obtained from Sulentic et al. (2006). 



ton humps are not accompanied by corresponding detections 
of Fe Ka lines. The results of the joint-fitting process also 
portray a mixed picture for the average reflection spectrum of 
luminous, high-accretion rate RQ AGNs as a class. The rel- 
atively weak constraints on the mean reflection properties of 
this class do not allow us to either confirm or rule out the exis- 
tence of an 'X-ray Baldwin Effect', regardless of whether the 
X-ray luminosity or the accretion rate drives the anticorrela- 
tion with EW(Fe Ka) (e.g., Zhou & Wang 2005; Bianchi et 
al. 2007). 

3.1.4. Soft Excesses 

By extending the X-ray spectral fitting to the entire EPIC 
energy range (0.2-10keV), we checked whether any of our 
nine high-redshift sources with XMM-Newton spectra (in- 
cluding the five high-redshift sources from S06) shows evi- 
dence for excess soft-X-ray emission, frequently observed in 
lower-redshift AGNs (e.g., Comastri et al. 1992; Reeves & 
Turner 2000; Piconcelli et al. 2005). While the physical na- 
ture of the soft excess is uncertain its presence is more pro- 
nounced among local high-accretion rate AGNs, i.e., narrow- 
line Seyfert 1 (NLS1) galaxies (e.g., Vaughan et al. 1999a; 
Boiler et al. 2002; Czerny et al. 2003; Vignali et al. 2004); 
hence it is of interest to search for its existence in our lumi- 
nous high-accretion rate sources at high redshift. We extrapo- 
lated the best-fit Galactic-absorbed power-law model obtained 
for rest-frame energies ^ 2 keV (see § l2.2l and Tabled to the 
0.2-10 keV observed-frame energy range and no signature of 
soft excess emission was detected in any of the sources (i.e., 
no systematic residuals were present for the extrapolated fits). 

This result is not unexpected given the relatively high red- 
shifts of our sources and the low-energy cutoff (0.2 keV) 
of the EPIC cameras. For example, Porquet et al. (2004) 
have found that the effective temperatures of the soft excess- 
components in a sample of moderate-luminosity sources at 
z < 0.5 are of the order of ~ 0.2 — 0.3 ke V. For our ten sources 
with z ~ 2, such temperatures shift to observed-frame values 
of ^0.1 keV. Furthermore, since the soft-excess component 
in most type 1 AGNs is typically observed below a rest-frame 
energy of ~ 0.7 keV (except for some high L/Lgdd NLSls, 
in which this component may extend up to 1.5 keV; e.g., 
Vaughan et al. 1999b), the detection of such a component with 
XMM-Newton is challenging even for sources with z ~ 1. 
Nevertheless, we have placed upper limits on the strength of 
a potential soft-excess component in two of our ten luminous 



sources that have the lowest redshifts, namely PG 1630+377 
and PG 1634+706 (with z = 1.476 and z = 1.334, respec- 
tively). In these two cases, the extension of the energy range 
to 0.2 keV has increased the number of photons in the fit 
by ~ 50%. The constraint on the data-to-model ratio in the 
2 keV rest-frame band, as a result of the extrapolation of 
the best-fit slope to the entire EPIC energy range, is ^,2.0 and 
£ 1.5, respectively. Finally, we note that the non-detection 
of soft-excess emission in PG 1247+267, PG 1630+377, and 
PG 1634+706 is consistent with previous studies of their 
XMM-Newton spectra (Page et al. 2004b; Piconcelli et al. 
2005). 

3.1.5. X-ray Variability 

Two of our sources, SBS 1425+606 and 
LBQS 0109+0213, were previously detected by the ROSAT 
PSPC and HRI, respectively. For the first of these, we mea- 
sured an unabsorbed flux of 4.7 x 10~ 14 ergs cm~ 2 s _1 in the 
0.5-2.0 keV observed-frame band from our XMM-Newton 
data, which is somewhat higher than, but consistent (within 
the errors) with the ROSAT flux (see also Just et al. 2007). 
Similarly, we measured for the second source an unab- 
sorbed flux of 2.0 x lO-^ergscm^s- 1 in the 0.5-2.0keV 
observed-frame band, which is a factor of ~ 3 lower than 
the value calculated from the ROSAT count rate. In both 
cases, we assumed that the photon index we measured from 
the XMM-Newton observations at rest-frame energies of 
^2keV also extends to the ROSAT bandpass (i.e., down 
to energies of 0. 1 ke V in the observed frame), based on th e 
lack of detectable soft excesses in these sources (§ 13.1.41 ). 
Although highly luminous RQ AGNs such as those presented 
here are not expected to exhibit pronounced X-ray variations 
even on timescales of several years (e.g., Lawrence & Pa- 
padakis 1993), the case of LBQS 0109+0213 is not unique, 
since several other luminous sources at high redshifts have 
displayed X-ray variations with similar amplitudes (e.g., 
Paolillo et al. 2004; Shemmer et al. 2005). 

We also searched for rapid (on timescales of ~1 hr in the 
rest frame) X-ray variations in the XMM-Newton and Chan- 
dra data of our five new sources by applying Kolmogorov- 
Smirnov tests to the lists of photon arrival times from the 
event files, but no significant variations were detected. To- 
gether with the results of S06, this suggests that high accretion 
rates do not necessarily lead to faster and higher-amplitude 
X-ray flux variations, as has been previously expected based 
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TABLE 3 

Compton Reflection and Iron Emission 



EW(Fe Reef 



AGN 




(eV) 


R b 


LBQS 0109+0213 


2.349 


<252 


3 - 3 -3.0 


2QZ 1023805.8-274337 


2.471 


<375 


< 1.4 


PG1247+267 C 


2.038 


< 188 


9 q+38.8 


Q 1318—113 


2.306 


<52 


1 4+3.0 


Q 1346-036 


2.370 


<298 


<7.2 


SBS 1425+606 


3.202 


< 613 


< 1.9 


PG 1630+377 


1.476 


458 +384 


<6.7 


PG 1634+706 c 


1.334 


<49 




HE 2217-2818 


2.414 


<233 


< 1.6 


NOTE. — Best-fit parameters of fitting each spectrum at the £ 2 keV 
rest-frame energy range with a model consisting of a Galactic-absorbed 
power-law, a Compton-reflection component, and a neutral Fe Kct emis- 



sion line. Errors represent 90% confidence limits, taking one parameter 
of interest (A^ 2 =2.71). 

a Rest-frame equivalent width of a neutral Fe Ka emission line at rest- 
frame E — 6.4 keV and a fixed width of o" = 0. 1 keV. 
k Relative Compton-reflection parameter; see text for details. 
c F-test indicates that the PEXRAV+ZGAUSS model results in an im- 
proved fit with respect to a power-law model. 

on X-ray variability properties of some NLSls, hence spoil- 
ing the analogy between NLSls and luminous, high-redshift 
AGNs (e.g., Grupe et al. 2006). On the other hand, our results 
are consistent with the idea that more massive BHs lead to 
longer timescales and smaller amplitudes of X-ray flux varia- 
tions (e.g., O'Neill et al. 2005). 

3.2. The Hard-X-ray Power-Law Photon Index 
3.2.1. Breaking the FWHMfHfi )-L/L Edd Degeneracy 

In Fig.g] we plot T versus FWHM(Hj3) and L/L Ed d for 
our core sample of 35 AGNs. Twenty five of these sources 
are Palomar Green (PG) quasars (Schmidt & Green 1983) 
at z < 0.5 with 44 < log[vL v (5100 A)] < 46, and ten are at 
z = 1.3 - 3.2 with 46 < log[vL v (5100 A)] < 48. In S06, 
it was shown that both FWHM(Hj3) and L/Lgdd are sig- 
nificantly correlated with F when only the 25 moderate- 
luminosity (z < 0.5) sources are considered (in agreement 
with Porquet et al. 2004; Piconcelli et al. 2005; see also 
Table [5}. This FWHM(Hj3)-L/LEdd degeneracy emerges as 
a consequence of considering only sources with a relatively 
narrow luminosity range, i.e., 44 < log[vL v (5100 A)] < 46 
(see also § 13.1.11 and S06). When the ten highly luminous 
AGNs at z = 1.3 — 3.2 are introduced to the analysis, both 
FWHM(H/3) and L/Lgdd remain significantly correlated with 
r (with > 99% confidence); however, the significance (in 
terms of chance probability) of the r-FWHM(H/3) correla- 
tion drops considerably while the significance of the F-L/LEdd 
correlation increases (Table [5}. The results of the F-L/LEdd 
correlation do not significantly change when a constant bolo- 
metric correction factor of 7 (which is the average correction 
factor for the luminosity range of our sample; see ^ 13. 1.1b is 
used to derive L/LEdd values; in this case, the significance of 
the correlation increases with the chance probability dropping 
from p = 1.6 x 10~ 3 to p = 4 x 10~ 4 . As suggested in S06, 
the extension of the luminosity range of the sample to high 
luminosities has allowed us to break the degeneracy between 
the r-FWHM(Hj3) and r-L/L E dd correlations and show that 
L/LE d d drives the correlation with F. 

We have also repeated the Mann- Whitney (MW) 
nonparametric rank test performed in S06 on all 
sources with FWHM(H/3) values that lie in the range 



TABLE 4 

COMPTON REFLECTION AND IRON EMISSION - JOINT FITTING 



Run 
No. 


No. of 
Sources 


No. of 
photons 




EW(Fe Ka f 
(eV) 


R 


I b 

n 

IV C 


9 
7 
5 
4 


39,173 
15,307 
9,814 
7,112 


2.22 
2.37 
2.39 
2.61 


<27 

< 114 

150+ 133 
L3V -109 

<300 


2 2 +L1 
^ -0.8 

"■'-0.6 

< 1.4 

< 1.1 



NOTE. — Best-fit parameters of joint fitting the spectra at the £ 2 keV rest- 
frame energy range with a model consisting of a Galactic-absorbed power- 
law, a Compton-reflection component, and a neutral Fe Ka emission line. 
Eirors represent 90% confidence limits, taking one parameter of interest 
(A^ 2 = 2.71). 

a The EW(Fe Ka) corresponds to a rest-frame at the given (z). 
b All sources from TablefJ] 

c All sources from TablefJ] excluding PG 1247 +267 and PG 1634+706. 

d All sources from Table fJJ excluding PG 1247+267, PG 1634+706, 

LBQS 0109+0213, and Q 1318-1 18. 

c All sources from Table fj] excluding PG 1247+267, PG 1630+377, 
PG 1634+706, LBQS 0109+0213, and Q 1318-118. 

3400<FWHM(Hj3)<8000km s" 1 [which is the FWHM(Hj3) 
interval of the ten luminous sources in our sample]. We find a 
significant deviation (with > 99.5% confidence) between the 
r values of the ten luminous sources and the nine moderate- 
luminosity (z < 0.5) sources in this range. In contrast, when 
the MW test is performed in the 0.1 L/LEdd ~ 0.5 range, 
the r values of the two groups of AGNs are not significantly 
different (as can be seen clearly from Fig. ]4}p). These results 
reinforce the S06 argument that F depends primarily on 
L/L E dd- 

We note that in the case of two of our luminous sources, 
LBQS 0109+0213 and SBS 1425+606, the F values are 
consistent with the expected values of both FWHM(Hj3) 
and L/LEdd, and hence they do not assist in removing the 
FWHM(Hj3)-L/L E dd degeneracy. The first of these sources 
has a remarkably flat hard-X-ray spectral slope (with F = 
1 .23; see § l2.2l and Table©, which was not expected given its 
original L/Lsdd determination from Shemmer et al. (2004; see 
their Table 2). Owing to the new and improved spectroscopic 
measurements of Netzer et al. (2007) for LBQS 0109+0213 
(SBS 1425+606), FWHM(Hj3) has increased from 5781 km 
s -1 to 7959km s _1 (3144km s -1 to 4964km s~ ! ), and hence 
L/L E dd has decreased from 0.2 to 0.1 (0.9 to 0.4). For the 
other Shemmer et al. (2004) sources used in this work, the 
Netzer et al. (2007) FWHM(H/3) measurements are consis- 
tent with those of Shemmer et al. (2004). The optical data 
presented in Tableware based on the new Netzer et al. (2007) 
measurements. 

Our result that F depends more strongly on L/L^dd, rather 
than on FWHM(H/3), may also serve as a partial vindica- 
tion for the use of the reverberation-mapping based method 
to determine L/Lsdd (see § 13.1. 11 1. From a pure mathemati- 
cal point of view, in the tests described above, we essentially 
multiplied one observational parameter, FWHM(Hj3), with 
another, vL v (5100 A), and yet the dependence of F on the 
product of the two parameters (i.e., L/LEdd) turned out to be 
stronger than the dependence on the single parameter. If the 
product of the two parameters had no physical significance, 
then there would be no apparent reason to expect this result. 

3.2.2. The T-L/Ledd Correlation 

We use the significant correlation we find between F and 
L/LEdd to derive a linear relationship between the two param- 
eters employing a variety of statistical methods. Hereafter, 
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FIG. 4. — The hard-X-ray photon index vs. FWHM(H/3) (a) and L/L Ed d (b). Circles mark sources at z < 0.5 with vL v (5100A)< 10 46 ergs s 1 from the 
Piconcelli et al. (2005) sample. Squares mark luminous sources at z = 1.3 — 3.2 with vL v (5100A)> 10 46 ergss -1 . Error bars on Y are shown at the 90% 
confidence level. 



TABLE 5 
Correlation Coefficients and 
Significance 



Correlation 


rs 


P 






N 


r-FWHM(HjS) 


-0.61 


1.2 x 


lo- 


3 


25 


r-FWHM(HjS) 


-0.44 


8.4 x 


ur 


-3 


35 


r-L/L Edd 


0.60 


1.5 x 


10" 


-3 


25 


T-L/Ledd 


0.55 


6.0 x 


10" 


-4 


35 



NOTE. — The last three columns represent the 
Spearman-rank correlation coefficient, chance prob- 
ability, and number of sources for each correlation, 
respectively. 

we take the error bars on F at the 1 a level. A standard % 2 
minimization method weighted by the errors on Y yields the 
following relation (with 1 a errors): 

r= (0.31 ±0.01)log(L/L Edd ) + (2.11 ±0.01). (1) 

The best-fit coefficients of this relation are consistent with 
those obtained by Wang et al. (2004) and Kelly (2007) 
who find a similar correlation in low-moderate luminosity 
sources. However, the % 2 value obtained by this minimization 
(% 2 / v = 980/33) suggests that either the data are not well 
represented by a linear model, the error bars on F are under- 
estimated, or there is additional intrinsic scatter in the data. 
Following the methods outlined in Tremaine et al. (2002) 
and K05, and by assuming a % 2 /v = 33/33 (i.e., a reduced 
% 2 = 1), we obtain an estimate of the additional potential scat- 
ter of AT ~ 0. 1 x r in the dependent parameter; this scatter is 
larger than the typical measurement errors on Y. 

The observed intrinsic scatter may be induced by the uncer- 
tainty in the BH-mass estimate, and/or additional unknown 
physical properties, such as the optical depth in the corona, 
orientation, and BH spin, that vary from source to source. An 
additional potential source for this scatter may be attributed to 
variability; for example, variability may seem to have a sig- 
nificant contribution to the scatter since the X-ray data and 
the rest-frame optical data (used for obtaining L/LEdd) are 
not contemporaneous. However, as discussed in S06, vari- 
ability (i.e., changes in both Y and L/LEdd) is n °t expected 
to dominate the observed scatter in this correlation (see also 
K05; Steffen et al. 2006; Wilhite et al. 2007); its effects are 



expected to be even less significant in this case, since our 
sources are mostly luminous and thus have smaller ampli- 
tudes of X-ray and optical flux variations (e.g., Lawrence & 
Papadakis 1993; Kaspi et al. 2007). 7 

To test whether Y can serve as an accretion-rate indicator we 
switched the roles of Y and L/LEdd, allowing the first to serve 
as the independent variable. To account for the scatter in the 
r-L/LEdd correlation, we performed a linear-regression anal- 
ysis using the Bivariate Correlated Errors and Scatter method 
(BCES; Akritas & Bershady 1996) on the data. For this pur- 
pose, we assumed that typical uncertainties on the determina- 
tion of Mbh are a factor of ~ 2 (e.g., Kaspi et al. 2000), and 
therefore assigned homoscedastic (i.e., uniform variance), la 
errors of 0.3 dex on log(L/LE dd ). The best-fit linear relation 
using the BCES bisector result (with 1 a errors) is 



log(L/L Edd ) = (0.9±0.3)r-(2.4±0.6). 



(2) 



A linear regression based on y} minimization using the FI- 
TEXY routine (Press et al. 1992) gives consistent results to 
those obtained with BCES with x /v = 34.24/33. In addi- 
tion, we also performed a linear-regression analysis on the 
data using the maximum-likelihood estimate (MLE) of Kelly 
(2007). The results of the different linear fits are presented in 
Tabled and are shown in Fig. [5] In all three linear-regression 
methods outlined above, the average scatter on the predicted 
value of log(L/LEdd) is ~ 0.35 dex. This is only slightly 
higher than the typical uncertainty, ~ 0.3 dex, associated with 
log(L/LEdd) determinations using the reverberation-mapping 
based Mbh determinations. 

We have also computed L/L^dd values using the Bentz et 
a l. (200 7) slope for the BELR size-luminosity relation (see 
§ 13.1. 11 1. This flatter slope (a = 0.54) results in M B n values 
that are smaller by a factor of ~ 0.7, on average, and hence 
L/LEdd values that are larger, on average, by the same factor. 
This has no noticeable effect on the strength of the r-L/L£ dd 

7 We note that some low-luminosity AGNs are known to exhibit large 
X-ray flux variations (with amplitudes of factors of ~ 10 or more) and corre- 
sponding X-ray spectral slope variations (with Ar ~ 0.4) that are often due to 
changing spectral contributions from reflection and/or absorption (e.g., Tay- 
lor et al. 2003; Krongold et al. 2007; Grupe et al. 2008). Such changes have 
not been observed for moderate-high luminosity sources such as those under 
study in this work. 
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FIG. 5. — Left: the L/Lem-T correlation (top) and residuals from the BCES fit (bottom). Symbols are identical to those in Fig.[4] Errors bars are at the lo" 
level. The solid line, dashed line, and dotted line mark the best-fit linear relation from the BCES, FITEXY, and MLE methods, respectively, outlined in the text 
(note that the three lines are almost overlapping; see Table[6). Right: distribution of the logfL/Ljidd) residuals from the BCES fit for the entire sample. 



correlation (as expected for a nonparametric correlation) or 
on its slope; only the intercept (i.e., a values in Tabled is af- 
fected, and it increases slightly within the current uncertainty 
on its value. We conclude that the F -L/LEdd relation is not 
significantly affected by small deviations from the K05 slope. 

Based on our high-quality sample of 35 sources, we con- 
clude that the hard-X-ray power-law spectral slope can pre- 
dict the value of the normalized accretion rate in RQ AGNs, 
across four orders of magnitude in AGN luminosity, with an 
acceptable uncertainty level of a factor of ,$,3. This may of- 
fer a useful new tool to probe the history of BH growth, based 
almost exclusively on the availability of high-quality hard- 
X-r ay s pectra of RQ AGNs; we discuss this in more detail 
in § l3.4l below. Although not unexpected, and previously pre- 
dicted by several studies (see §[T), our result is by far the most 
reliable indication that the shape of the hard-X-ray power-law 
spectrum is largely controlled by L/LEdd- As mentioned in 
§[T] a possible explanation is that the corona acts as a 'ther- 
mostat' by cooling more efficiently when the disk emission 
increases, manifested by a steepening of the hard-X-ray spec- 
trum. 

In this context, it is also interesting to note that our re- 
sult may readily explain the narrow ranges observed for the 
values of both F and L/LEdd in optically selected moderate- 
high luminosity, type 1 RQ AGNs. While L/LEdd is rela- 
tively narrowly distributed around a value of ^0.3 (0.28 in 
our sample of 35 sources), with a typical dispersion of a fac- 
tor of ~5 around that value (e.g., McLure & Dunlop 2004; 
Kollmeier et al. 2006; Netzer et al. 2007; Shen et al. 2008), 
r values for such sources typically lie in the narrow range 
of ~ 1.5-2.5 (1.2-2.6 for our sample of 35 sources; see also 
Vignali et al. 2005). These ranges in L/LEdd an d T may be dif- 
ferent for lower luminosity sources. In particular, for sources 
with vL v (5100A) ;$ 10 42 ergss~', a possible difference may, 
in part, be attributed to a difference in their accretion mode as 
compared with the moderate-high luminosity sources studied 
in this work. Finally, we point out that a possible dependence 



of the bolometric correction for Z.2-i0keV on L/Ledd, recently 
reported by Vasudevan & Fabian (2007), may, in part, be re- 
flected by our results. 

3.3. What Determines the Optical-X-ray SED? 

The a ox parameter is known to have a strong anticor- 
relation with optical-UV luminosity at 2500 A [hereafter 
L V (2500A)], and it shows no significant dependence on 
redshift (e.g., Vignali et al. 2003; Strateva et al. 2005; 
Steffen et al. 2006; Just et al. 2007; but see also Kelly 
et al. 2007). This is equivalent to a non-linear rela- 
tion between X-ray and optical-UV luminosity of the form 
L v (2keV) oc L v (2500A) a , where a < 1 (see Eq. 7 of Just et 
al. 2007 who find a = 0.709 ± 0.010). This relation still lacks 
a sound physical interpretation and, in particular, it is not clear 
what mechanism controls the proportion of reprocessed hard- 
X-ray emission from the corona with respect to the UV emis- 
sion from the disk. 

In order to test whether Oo X depends on L/Lgdd^ we selected 
8 1 sources from the Steffen et al. (2006) sample with available 
FWHM(H/3) measurements (all these have z ^ 0.75 and they 
i nclude the 25 sources with z < 0.5 from our core sample; 
§ 13.2. Il l, and complemented this sample with the ten luminous 
sources at z = 1.3 — 3.2 from S06 and this work. For about 
half of the Steffen et al. (2006) sub-sample, FWHM(Hj3) 
values were obtained from Boroson & Green (1992), and 
for the rest we used FWHM(Hj3) measurements from Net- 
zer & Trakhtenbrot (2007), kindly provided by B. Trakhten- 
brot, 2007, private communication. For each source we deter- 
min ed th e value of L/Lgdd following the pro cedure outlined 
in § 13.11 Following the arguments made in § 13.2.21 the non- 
simultaneous X-ray and optical observations are not expected 
to affect our analysis considerably (see also § 3 of Steffen et 
al. 2006). 

We find that a ox is significantly correlated with both 
L V (2500A) and L/Lgdd; however, the correlation with lumi- 
nosity is considerably tighter (see Fig. [6). Even when a con- 
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TABLE 6 

Linear Regression Coefficients for the L/L Ed d-r correlation 



BCES bisector FITEXY MLE a 

P a p a p a 

0.93±0.31 -2.44±0.63 0.94±0.21 -2.46±0.42 0.97 ±0.26 -2.50±0.50 
NOTE. — The slope and intercept for each method are represented by /J and a, respectively. Errors are 

at 1 o" confidence levels. 

a See Kelly (2007) for more details. 



stant bolometric correction facto r of 7 is used to determine 
L/Lsdd values (as performed in § 13. 2. Il l, the OC OX -L V (2500 A) 
correlation is considerably tighter than the a ox -L/Ledd cor- 
relation, although the significance of the latter correlation 
has increased (the chance probability has dropped from 
p = 1.3 X 10~ 3 to p = 1.3 x 10~ 4 ). By investigating this fur- 
ther, we find that the a x-L/Ledd correlation disappears when 
the Oo X -L v (2500 A) correlation is taken into account; i.e., the 
difference between the observed a ox and the expected a ox 
[based on the L V (2500A) values of the sources and the most 
recent C£ ox -Lv (2500 A) relation from Just et al. 2007] is not 
correlated with L/Ledd- In another test, we divided the data 
into sources that have L/Ledd values either above or below the 
median value of L/Ledd= 0.22. In both sub-samples we de- 
tect a strong Oo X -L v (2500 A) correlation (and no Oo X -L/Ledd 
correlation). These tests suggest that the a x-L v (2500 A) cor- 
relation is persistent and that Oo X does not depend primarily 
on L/Ledd- In this scenario, the (relatively weak) Oox-L/Ledd 
correlation may be a consequence of the inherent dependence 
of L/Ledd on L V (2500A), and the relative weakness of this 
correlation compared with the strong a x-L v (2500 A) corre- 
lation may be due to the additional uncertainty introduced by 
determining L/Ledd (see, for example, the typical error bars in 
Fig. [6]). Subsequently, if a O x-L v (2500 A) is the fundamental 
correlation, then the inclusion of the FWHM(Hj3) parameter 
to produce L/Ledd adds substantial scatter resulting in a rather 
weak a ox -L/LEdd correlation. It is also of interest to point 
out that we find no correlation between a ox and F (that de- 
pends on L/Ledd) m our core sample of 35 sources discussed 
in §E2J 

We also checked whether the 0Co X -L v (2500 A) correlation is 
induced by luminosity-dependent obscuration (e.g., Lawrence 
1991; Gaskell et al. 2004). In a first test, we found a strong 
correlation between Oo X and L v (5100 A) for our sample, with 
similar coefficients as in the a ox -L v (2500 A) correlation. We 
note that the fluxes at 5 100 A were obtained independently 
from the fluxes at 2500 A, i.e., no extrapolations have been 
made between the flux densities in these two wavelengths. 
In a second test, we redefined Oo X by changing the optical- 
UV continuum threshold from 2500 A to 5 100 A; the corre- 
lations between this modified a ox and both L V (5100A) and 
L V (2500A) for our sample returned similar correlation co- 
efficients as in the oc^-Ly (2500 A) case. We conclude that, 
within the limits of our sample, it is unlikely that the a ox - 
L V (2500A) correlation is induced by potential reddening ef- 
fects. 

The fact that L v (2500 A) appears as the underlying parame- 
ter determining the slope of the optical-X-ray SED instead of 
L/Lgdd is puzzling, since L/Ledd indicates the relative AGN 
power and significantly affects the SED of the disk. One pos- 
sibility is that our empirical expression for L/Ledd (see S06) 
is inaccurate or more complicated. This, however, is not sup- 



ported by the fact that F does depend on L/Ledd as defined 
here and not on vL v (5100A) (see §[T]and §1221] S06). We 
have also considered the possibility that the a 0X -L v (2500 A) 
correlation is induced by selection effects (e.g., Green et al. 
2006), but this is highly unlikely given the fact that the corre- 
lation spans over five orders of magnitude in luminosity (see, 
e.g., Just et al. 2007). Indeed, future correlations between 
a ox , L v (2500 A), and L/Ledd should inv olve L/Ledd determi- 
nations for X-ray-selected sources (see § l3.4l below). 

An alternative approach to interpreting the non-dependence 
of Oo X on L/Ledd stems from the recent results of Vasude- 
van & Fabian (2007). These authors confirm the strong de- 
pendence of Oo X on L V (2500A) in their AGN sample, but 
they do not find an 0C O x-L/Ledd correlation. Utilizing far-UV, 
as well as optical, near-UV, and X-ray spectra for a sample 
of AGN with known Mbh and L/Ledd. they argue that the 
reason for the non-dependence of a ox on L/Ledd lies in the 
choice of the optical-UV and X-ray continuum thresholds 
(i.e., 2500 A and 2keV) for the Ofex definition. While they 
find dramatic differences between the optical-X-ray SED of 
sources with different L/Ledd values, the ratio between the 
flux densities at 2500 A and 2 ke V remains almost constant as 
L/Ledd changes. The most pronounced differences between 
the optical-X-ray SEDs of different sources are expected to 
be concentrated in the m 0.01 — 0.1 keV spectral region (i.e., 
the 'big blue bump'), which cannot be traced effectively us- 
ing the current optical-UV and X-ray thresholds of a ox . In 
this context, the a ox parameter, given its current definition, 
while providing useful information on the Luv-Lx connec- 
tion, cannot be used as an accretion-rate indicator. By in- 
spection of Fig. 13 of Vasudevan & Fabian (2007), one may 
expect a strong dependence of oc ox on L/Ledd if, for instance, 
the optical-UV threshold is shifted from 2500 A to as close as 
possible to ~ 250 A, although this is observationally challeng- 
ing and may not be practical. When using modified definitions 
for a ox , there also remains the question of which correlation, 
O£ ox -L/Ledd or a ox -L\jy, is stronger than the other. In general, 
the choice of optical-UV and X-ray thresholds should likely 
depend upon the scientific question of interest. 



3.4. 77ie AGN Hard-X-ray Spectrum as a Probe for BH 
Growth 

In ij |3.2.2l we have shown that F may serve as an L/Ledd in- 
dicator in unabsorbed, moderate-high luminosity RQ AGNs 
with an acceptable mean uncertainty level on L/Ledd of a 
factor of 3. Therefore, in principle, given a high-quality 
X-ray spectrum in the ^ 2 keV rest-frame band of a RQ 
AGN, one can empirically estimate L/Ledd and Mbh for the 
source. Accurate measurements of F and Lx (in the 2-10keV 
rest-frame band, for example) may provide L/Ledd and Lt>d 
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FIG. 6. — Correlations between a ox and (a) L V (2500A) and (b) L/L^dd- Symbols are similar to those in Fig. [4] although the low-moderate luminosity sources 
were obtained from Steffen et al. (2006). The solid line marks the best-fit relationship between a ox and L V (2500A) from Just et al. (2007), and the cross marks 
the typical uncertainties on ce ox and L/L^dd- Spearman-rank correlation coefficients and chance probabilities are indicated in each panel. 



using our Eq.|2] and Eq. 21 of Marconi et al. (2004), 8 re- 
spectively, and thus provide an estimate of Mbh. We note 
that the Marconi et al. (2004) determination of Lbol (whether 
from L opt or Lx) relies, in part, on the very strong and 
non-linear Lx-iuv dependence (see § I3.31 >. The observed 
(luminosity-dependent) rms errors on that relation are given 
in Table 5 of Steffen et al. (2006); these represent typical 
deviations of up to ~ 0.15 from the mean 0Co X [i.e., a fac- 
tor of ~ 2.5 uncertainty on (/2kev//2500 a)] f° r our lumi- 
nosity range. This scatter is inherent in the L/Le^-F cor- 
relation as it reflects uncertainties in determining Lbol from 
vL v (5100A) and potential X-ray-optical variability. For ex- 
ample, when Mbh values for our core sample of 35 sources are 
recovered from L/Lgdd using the bolometric corrections for 
L2-iokeV, trie rat i° between t hese M bh values and the orig- 
inal Mbh values obtained in § 13.1. H is 1, on average, with a 
dispersion of a factor of 2. In an additional test, we used 
the r and L2-i0keV values of six nearby sources from our 
sample to recover their Mbh values and compare them with 
the most recent, direct (reverberation-mapping based) Mbh 
measurements from Peterson et al. (2004). The six sources, 
namely PG 0804+761, PG 0844+349, PG 0953+414, 
PG 1211 + 143, PG 1307+085, and PG 1613+658, comprise 
~ 40% of all the RQ AGNs with vL v (5 100A) > 10 44 ergs s" 1 
in the Kaspi et al. (2000) sample (i.e., the PG quasar sample) 
that have direct Mbh measurements. We find that the ratio be- 
tween the recovered X-ray-based masses and the correspond- 
ing Mbh measurements from Peterson et al. (2004) is 0.8, on 
average, with a dispersion of a factor of 1.9; we also find that 
the mass recovered for each individual source is consistent 
with the measured value, given the uncertainties from Eq. [2] 
and those from the reverberation-mapping measurements. 

High-quality X-ray spectra may thus be useful for estimat- 
ing the accretion rates and BH masses for moderate-high 
luminosity type 1, RQ AGNs and, in particular, for X-ray- 

8 Vasudevan & Fabian (2007) claim that the bolometric correction factor 
required to transform L2-10 keV mt0 A>ol depends on L/Z.Edd. while Marconi 
et al. (2004) have used the Vignali et al. (2003) Oo x -L V (2500A) correlation 
(which is not significantly different from the most recent correlation of this 
type given by Just et al. 2007) to derive bolometric corrections that depend on 
Z,(j i. Accurate measurements of T (that provides L/Z.Edd) an d will enable 
comparisons of L^„[ values obtained from these two methods. 



selected sources. This method may allow tracing the his- 
tory of BH growth in the Universe by utilizing large AGN 
datasets (see e.g., Brandt & Hasinger 2005) in which BH 
growth cannot be determined effectively using existing meth- 
ods [i.e., reverberation-mapping based methods (aka single- 
epoch methods), e.g., K05, and host-AGN type relations, e.g., 
Marconi & Hunt (2003)], either due to the faintness of the 
sources in the optic al-near-IR bands, or that the required 
spectroscopic features are either not accessible or cannot be 
modeled reliably. One advantage of the X-ray-based method 
is the ability to obtain L/L^dd an d Mbh by measuring broad 
spectroscopic properties (i.e., F) as opposed to detailed spec- 
troscopy for obtaining the width of an emission line (e.g., 
Hj3 or Mg II; this involves, for example, careful decontamina- 
tion of Fe II emission features from the UV-optical spectra). 
Moreover, measurements of the hard-X-ray power-law spec- 
tral slope are typically not limited to specific redshift ranges, 
such as those dictated by atmospheric transmission bands (and 
detector bandpass) for ground-based spectroscopy. These ad- 
vantages offer a way to obtaining many L/Lgdd an d Mbh 
estimates, economically, in contrast with more complicated 
(and redshift-restricted) spectroscopic measurements for in- 
dividual sources done with existing methods. Nevertheless, 
we caution that the X-ray measurements of F and Lx should 
be done carefully to account for potential complex absorption 
and Compton reflection. 

The X-ray-based method for estimating L/Lgdd an d Mbh 
may prove to be even more rewarding in cases where a source 
is either optically faint and/or the broad-emission lines are 
too weak to measure if, for example, these lines are over- 
whelmed by host-galaxy continuum (e.g., Moran et al. 2002; 
provided the source is not radio loud and that the redshift 
can be determined). In particular, this method may be the 
best way to determine L/Lgdd an d Mbh directly in obscured 
(i.e., optical type 2), moderate-high luminosity AGNs. Pro- 
vided the X-ray absorption column density is not too high 
(N R < 10 23 cm~ 2 ) and can be modeled accurately, the intrin- 
sic hard-X-ray power-law spectrum can, in principle, provide 
L/Ledd and Mbh a s outlined above. Examples of such sources 
with high-quality X-ray spectra are given in e.g., Civano et al. 
(2005), Mateos et al. (2005), and Mainieri et al. (2007). Since 
Compton-thin type 2 AGNs comprise a significant fraction of 
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the AGN population, tracing the growth of the supermassive 
BHs in their centers using the X-ray-based method is crucial 
for testing models of the evolution of BH growth and accre- 
tion luminosity in the Universe in an unbiased way. 

4. CONCLUSIONS 

We present X-ray spectroscopy for five highly luminous RQ 
AGNs at z = 1.3 - 3.2, with accurate FWHM(Hj3) measure- 
ments that allow determinations of their Mbh and L/Lem val- 
ues. Analysis of the X-ray spectra provided measurements of 
the hard-X-ray photon index in the rest-frame £ 2 keV band 
and a ox . We have combined these data with the S06 sample 
of 30 moderate-high luminosity sources with similar proper- 
ties, while doubling the number of highly luminous sources 
in their sample. Our main goal was to test the S06 claim that 
r can serve as an accretion-rate indicator in RQ AGNs. We 
have also tested whether any additional X-ray properties of 
our highly luminous sources depend on L/Ledd- Our main 
results are summarized as follows: 

1. Our new highly luminous sources with FWHM(Hj3) 
measurements have allowed us to break the degen- 
eracy between the dependence of F on FWHM(Hj3) 
and on L/Ledd, suggesting that the accretion rate 
largely determines the hard-X-ray spectral slope across 
four orders of magnitude in AGN luminosity (i.e., 
44 < log[vL v (5100 A)] < 48). 

2. We found a significant correlation between L/Lem an d 
r with a best-fit line of the form log(L/LEdd) = (0-9 ± 
0.3)r — (2.4 ± 0.6), and an acceptable uncertainty of a 
factor of £ 3 on a predicted value of L/Lem- 

3. Utilizing a sample of 91 sources from Steffen et al. 
(2006) and this work, we find that a ox depends strongly 
on optical-UV luminosity and only weakly on L/Leaa; 
the (weak) correlation with L/Lgdd is probably due to 
the strong L-L/Leu dependence. We discuss possible 
explanations for this result including the possibility that 
CCqx cannot be used as an accretion-rate indicator based 
on its current definition. 

4. We find a significant Compton-refiection feature in two 
of our sources, and the mean relative reflection for 



seven other sources is R = 0.9^Q g. By setting rather 
loose constraints on the strengths of Fe Ka emission 
lines in our highly luminous sources, we can neither 
confirm nor rule out a suggested anticorrelation be- 
tween EW(Fe Ka) and either luminosity or L/Leaa', 
the upper limit on the mean rest-frame EW(Fe Ka) 
for seven of these sources that do not show Compton- 
refiection features is 105 eV. 

5. We have not detected any signature of a soft-excess 
component in any of our highly luminous sources, in- 
cluding two sources at z ~ 1.4 where our rest-frame 
coverage extends to ~ 0.7 keV, suggesting that the soft 
excess does not depend strongly on the accretion rate. 

6. Although one of our highly luminous sources, 
LBQS 0109+0213, exhibits long-term (i.e., on 
timescales of years) X-ray variations, rapid X-ray vari- 
ations on timescales of ^lhr have not been detected 
in any of our highly luminous (and high-Men) sources, 
supporting the idea that X-ray variability timescale 
depends inversely on Mbh and does not depend on 
L/L E dd- 

The strong correlation between F and L/Lem may serve as 
a useful probe for tracing the history of BH growth in the uni- 
verse. It may provide L/Leaa and Mbh estimates for X-ray- 
selected AGNs, with the possibility of estimating these prop- 
erties for Compton-thin type 2 AGNs for the first time. 
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